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Abstract. - The critical fluctuations of superconductors are discussed in a fixed dimension 
scaling suited to describe the type II regime. The gauge dependence of the anomalous dimension 
of the scalar field is stablished exactly from the Ward-Takahashi identities. Its fixed point value 
gives the r\ critical exponent and it is shown that r\ is gauge independent, as expected on physical 
grounds. In the scaling considered, n is found to be zero at 1-loop order, while v » 0.63. This 
result is just the 1-loop values for the XY model obtained in the fixed dimension renormalization 
group approach. It is shown that this XY behavior holds at all orders. The result rj = ijxy 
should be contrasted with the negative values frequently reported in the literature. 



The high temperature superconductors have a larger critical region relative to ordinary 
superconductors. This fact allows in principle an experimental access to the critical region and 
measurements of critical exponents had been made in these materials, specially the YBCO 
ones jj], 0]. The experiments showed that critical fluctuations are very important in high 
temperature superconductors, a consequence of this fact being the non- mean- field values of 
the critical exponents. In the tested temperature region it was obtained that —0.03 < a < 
and v sa 0.67 while the amplitude ratio A + /A~ = 1.065 ± 0.01. These values are consistent 
with those found for He and correspond to a uncharged d = 3 XY universality class. There 
are some subleties with respect to this conclusion and the interpretation of the experimental 
data. For instance, the value of the v exponent has been obtained by a direct measurement 
of the London penetration depth, A, which was found to scale as A ~ t~ y , t being the reduced 
temperature and y = 0.33 ± 0.01 0. The value of v follows by assuming that the superfluid 
density scales with A as p s ~ A -2 and through the use of Josephson's relation p s ~ ||. It 
has been argued by Herbut and Tesanovic M that in a superconducting transition A should 
diverge with the same exponent as £, contradicting the relation p s ~ A~ 2 . In fact, they claim 
that the correct relation should p s ~ A -1 , when the critical fluctuations near the charged fixed 
point are taken into account. Thus, it seems to be very difficult to have an experimental access 
to the charged fixed point, even in the case of cuprates superconductors. The critical region 
probed in Refs. pi || corresponds to a crossover regime where the gauge field fluctuations 

Typeset using EURO-M45X 



2 



EUROPHYSICS LETTERS 



are unimportant. Theoretically, the nature of the charged fixed point has been elucidated in 
the early eighties through the study of lattice abelian gauge models using duality arguments 
[H. It has been shown that the normal-superconducting transition should be a second order 
phase transition, at least for type II superconductors. This means that it must exists an 
infrared stable charged fixed point, contradicting the weak first order transition scenario of 
the e-expansion [ [To| . More recent studies performed directly in continuous models |^, [?], ^) 
and recent numerical simulations in the lattice ||, gives further support to this view. 

However, some important theoretical questions concerning the scaling in the Ginzburg- 
Landau (GL) model are still open. For instance, everybody agrees that the critical behavior 
should be the XY one, but nobody is able to find the XY value for the 77 exponent from 
renormalization group calculations. More importantly, the value generally found in these 
calculations is negative. Negative values of 77 does not violate the scaling relations provided 
77 > — 1 and, in fact, slightly negative values had been reported and this bound is fulfilled 
@ 0, g ]l| § g|. However, the best XY estimate is slightly positive, 77 « 0.04. 
Therefore, the negative values found are possibly an artifact of perturbation theory. It is 
worth to mention, however, that respectable values for v consistent with the XY behavior are 
found in some RG calculations ^ |L §1. 

Nearly five years ago, Kiometzis and Schakel (l5j argued that negatives values of 77, though 
> — 1 , should be unphysical since in principle it violates unitarity in the corresponding quantum 
field theory Moreover, 77 is the fixed point value of the anomalous dimension 77^ of the order 
parameter field in the GL model, a quantity which is gauge dependent. Based on physical 
intuition only we may suspect that 774, should be gauge independent at the fixed point. However, 
there is no actual proof of this fact to date. For this reason we shall address this point in this 
paper by using a fixed dimension RG approach for T > T c . As will be made apparent soon, 
the critical point fixed dimension approach employed in references || is not well suited 
for the following analysis. By making use of the Ward-Takahashi (WT) identities, we shall 
establish the gauge dependence of 77^ at all orders in perturbation theory. This step follows 
from general field theoretical arguments |[6| . Next we shall show that at the critical point 77^ 
is in fact gauge independent. Finally, we perform a 1-loop calculation to obtain 77 — 0, a result 
consistent with the 1-loop result for the XY model. We show that higher order corrections 
will improve further this result to obtain 77 = tjxy- 

Our starting point is the bare action for the GL model, or Euclidean scalar QED in three 
dimensions, 



S = d A x 



+ Sgf, (1) 



where the zeroes denote bare quantities, Fq is a short for Fq U Fq U and £)° = <9 M + ie^A^. The 
S g f is the gauge fixing part and is given by 



Sgf = J d 3 x^-(d,Al) 2 . (2) 

We introduced a mass to the vector field in order to regularize the infrared divergences 
arising from some diagrams containing gauge fields propagators. This mass term breaks gauge 
invariance but as we shall see, gauge invariance is restored at the infrared stable fixed point. 

The renormalized action is defined by S' + SS where S' is the same as the bare action but 
with renormalized quantities while SS is the counterterm action. It is given by 
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with the renormalized fields defined by A^ = Z 

By adding sources terms for the corresponding fields, it is straightforward to derive the 
following WT identity: 



M 



-A i n 



-J{x) 



SJ{x) 



W(J IM ,j\j)=dMx), (4) 



where W = log Z , Z being the generating functional of correlation functions. When the sources 
are zero Z corresponds to the partition function. The W generates the connected correlation 
functions. The Legendre transform of W is performed as usual and gives the functional 
T(ip', <p } a ) which is the generator of the 1-particle irreducible functions. It satisfies a WT 
identity which is the Legendre transform of (0): 



-A - M' 
a 
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(5) 



The WT identity given by Eq.(^J) gives important informations about the counterterms. 
For example, it implies that non-gauge invariant terms in the renormalized action are not 
renormalized and consenquently the corresponding counterterms are zero. This implies M 2 = 
ZaMq and a = Z^cto- Gauge invariance also implies e 2 = Za€q- Let us define the following 



dimensionless gauge couplings, e 2 
equations: 



e /m and v = m/M. We have the following exact flow 
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(6) 
(7) 
(8) 
(9) 



where we have introduced the RG function tja which is the anomalous dimension of the gauge 
field. It is defined by 



VA=m— log Z a- (10) 
dm 

An immediate consequence of (B) is that at a charged fixed point we must have 77^4 = 1. This 
very simple observation have important implications concerning the scaling of the magnetic 
field penetration depth, as first observed by Herbut and Tesanovic |J (see also ref. ]l7|). 
It follows also from the above equations that at the charged fixed points the corresponding 
fixed point value of a is a* — 0, that is, the Landau gauge. Note that the situation here is 
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somewhat different from that one encountered in particle physics where d = 4. In fact, in that 
case the beta function for e 2 is t\a e 2 and, therefore, the charged fixed point would correspond 
to t\a = 0. Since Eq.(|^) remains the same for d = 4, we obtain that the fixed point value of a 
is arbitrary in this case. 

From Eq.(Q) we obtain that the charged fixed point value of M 2 is zero. Therefore, near the 
superconducting fixed point the effective action flows to a configuration with massless gauge 
fields in the Landau gauge. This explains why RG calculations performed in the Landau gauge 
gives good results. We can say, therefore, that it is legitimate to compute critical exponents 
in the Landau gauge even knowing that some RG functions like r/^ = mdlogZ^/dm are gauge 
dependent. This is in contrast to 77^4 which is gauge independent if a minimal subtraction 
scheme is used. 

In order to check the consistency of this argument it remains to show that 77^ is in fact 
well behaved with respect to the gauge dependence as approaching the critical point. We can 
obtain the exact gauge dependence of 77^ from the WT identity Eq.(|J). Indeed, we can use 
(||) to relate scalar 2-point correlation functions with a (d^A^) 2 insertion to the 2-point scalar 
correlation functions. This is obtained by applying twice the WT identity (||) to the 2-point 
correlation functions. We obtain in this way the exact equation 



W, 



(2 



d 3 k 



(2tt) 3 (k 2 +aM 2 ) 2 
out of which we get the exact gauge dependence of 77^: 



W {2 \p + k) - W {2 \p) 



(11) 



= rj A e 2 v 
da Siry/a' 

with t]a being gauge independent. As a charged fixed point is approached, t]a — > 1, v — ► 0, 
e — > e* and a — ► 0. In order to have a consistent 77 exponent it is necessary to show that 
v/y/a — » as we approach the critical point. This is in fact the case since from Eqs.dTl) 
and (§) we obtain that v/y/a scales like m = £ 1 near the critical point. Therefore, we 
have drjfi/da — at the critical point. The same result is obtained if we consider a more 
general model including a Chern-Simons term whose critical behavior has been studied recently 
|§ |l|, |||. Although in this case the photon has a massive propagator without breaking of 
gauge invariance, we must still keep the mass M above in order to recover consistently the 
non-topological model in the zero Chern-Simons mass limit. Note that the critical point 
approach B|, though it regularizes infrared divergent graphs, is not appropriated to discuss 
in a complete way the gauge dependence of 77^. In fact, at the critical point m — M — and 
we have infrared divergences in (|IT|). 

Let us perform now a sample calculation up to 1-loop order in the Landau gauge. For this 
end we write <j) — (4>i + ifo) / V% and use the the renormalization conditions for the irreducible 
vertex functions: T^(0) = to 2 , r$ = 3M 2 , dT^(0)/dp 2 = 1, dT^{Q)^/dp 2 = 2 and 
liiii(O) = 3m. The corresponding anomalous dimensions are given by 



flA - ^ (13) 



2 

3^r(l + ?j) 

The flow of the coupling u = u/m is given up to 1-loop order by 



^ = ~— H , 32 " ( 14 ) 
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m ^L = {2^-l)u+^-e + ^- e\ (15) 

Figures 1 and 2 show the flow diagram respectively in the (u, /) (/ = e 2 ) and the (/, v) planes. 
Fig.l corresponds to a section v — 0.001 of the critical manifold. Note that for a small but 
nonzero v we have two charged fixed points, corresponding respectively to the tricritical and 
superconducting fixed points Jfi, 14, §|. It is useful to compare the above calculation with other 



fixed dimension approaches, for instance, the one employed by Herbut and Tesanovic and de 
Calan et al. In the approach of Refs. || || the charged fixed points are obtained in a 

critical point calculation through the introduction of a constant parameter c, giving the ratio 
between two different momentum scales of the problem, namely, the momentum determining 
the running of e 2 and the one determining the running of ii. The parameter c can be adjusted 
in order to generate charged fixed points. The arbitrariness of c is removed by fixing it from 
a known numerical value of the Ginzburg constant, k, at the tricritical fixed point. In our 
case, it is v that plays the role of c, v representing the ratio between the two existing scales 
in our problem, namely, m and M. An important difference between the present approach 
and the one of Refs. || || is that we do not need to fix numerically v since it flows naturally 
to a fixed point value. Note that the infrared stable fixed point in the flow diagram of Fig. 2 
is charged. The problem with the e-expansion is that only one scale is considered. In such 
a model we have naturally two scales, which are of course related. For instance, the most 
natural scales in the broken symmetry phase are A and £, whose ratio gives an important 
physical parameter, the Ginzburg constant k. However, the e-expansion has the advantage of 
being a controlled approximation, in the sense that we have a well defined small parameter. 
In the fixed dimension approach, a good choice of expansion parameter is 1/N, N being the 
number of order parameter components. The main drawback in this case is that it is not easy 
to extrapolate the value of N to the physical case N = 2. Thus, our approximation, though 
uncontrolled (just like those in Refs. |g, |[ [li], enables us to get sensible physical results. 
Anyway, it is possible in principle to use controlled approximations like the e-expansion to 
obtain results consistent with the existence of an infrared stable charged fixed point. This 
can be accomplished by considering explicitly the two scales of the model. This problem is 
treated more appropriately in the broken symmetry regime and will be the subject of a future 
publication. 

Since v = at the charged fixed point, we obtain exactly the critical behavior of the XY 
model. For instance, we find r\ = and v 0.63, the 1-loop values for XY model in the 
fixed dimension approach. Note that r\ is not negative in the present scaling. Higher order 
corrections behave in the same way, that is, all the powers of e 2 are suppressed at the charged 
fixed point because they are multiplied by some function of v which goes to zero at the charged 
fixed point. Note that these functions of v will never be divergent as v — > 0. This follows by 
simple dimensional analysis performed in the graphs containing gauge and scalar fields lines. 
Thus, if we compute higher order corrections we expect to approach asymptotically the best 
XY values estimates for r\ and v. 

In summary, we stablished exactly the gauge dependence of the scalar field anomalous 
dimension and showed that calculations are legitimate if performed in the Landau gauge. In 
the fixed dimension scaling considered, explicit calculations show that the GL model lies in 
fact in the XY model universality class. 
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Fig. 1. - Flow diagram in the (tt, /)-plane. 
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Fig. 2. - Flow diagram in the (/, u)-plane. The infrared stable fixed point is charged 
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